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Abstract 

We study the destabilization mechanism in a unidirectional ring of identi¬ 
cal oscillators, perturbed by the introduction of a long-range connection. It is 
known that for a homogeneous, unidirectional ring of identical Stuart-Landau 
oscillators the trivial equilibrium undergoes a sequence of Hopf bifurcations 
eventually leading to the coexistence of multiple stable periodic states resem¬ 
bling the Eckhaus scenario. We show that this destabilization scenario persists 
under small non-local perturbations. In this case, the Eckhaus line is modu¬ 
lated according to certain resonance conditions. In the case when the shortcut 
is strong, we show that the coexisting periodic solutions split up into two 
groups. The first group consists of orbits which are unstable for all parameter 
values, while the other one shows the classical Eckhaus behavior. 


1 Introduction 

The theory of dynamics on networks is a rapidly growing field of research. Motivated 
by applications from nenroscience, biological systems, and social networks, scientists 
from varions areas have investigated the connection between network strnctnre and 
dynamics in the last few decades. Since then, the importance of special topolo¬ 
gies like scale-free or small world networks was discovered and dynamics on these 
strnctnres were stndied extensively [1]. 

In this article, we stndy a special class of network topologies, namely nnidirec- 
tional rings which are pertnrbed by the insertion of a single additional non-local 
link [Fig. 1]. One reason for the interest in ring strnctnres is that they emerge in 
many natnral systems [2-4]. Also, rings can be seen as motifs of larger and more 
complex networks [5]. A lot of research has been done for bidirectional rings in the 
last few years [6-10]. In contrary, less is known abont dynamics in nnidirectional 
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rings [11-16], although these structures play an important role in various applica¬ 
tions [4, 17-21], As a simple, paradigmatic model, we consider unidirectional rings 
of N identical Stuart-Landau oscillators with an additional shortcut from node i 
to node N, to which we refer as a perturbation of the homogeneous system. The 
dynamics on the perturbed ring are described by the following equations: 


Zj (t) = (/i - \zj {t)f) Zj + Zj+i (t), j = 1, ..., - 1, 

Zn (t) = (/i - \zn (t) 1^) Zn + (t) + SZi (t) , 


( 1 ) 


where ^ = a + i(3, (3 > 0, s > 0 is the shortcut strength, and Zj (t) G C. For an 
unperturbed ring (s = 0), the destabilization mechanism is described in [14], In 
particular, when the parameter a = (/i) in (1) is increased, the zero equilibrium 


( 2 ) 


Zi = ... = Zn = 0 


loses stability and undergoes a sequence of Hopf bifurcations. The first half of 
the emerging branches of periodic solutions stabilizes at appropriate values of a 
(see chapter 2.2). Remarkably, this phenomenon resembles the Eckhaus scenario in 
spatially extended diffusive systems [22, 23]. The aim of this article is to investigate 
the transformation of this scenario under non-local perturbations which destroy the 
rotational symmetry of the system. We investigate two different asymptotic cases 
of small and large perturbation size s. For small s, the Eckhaus scenario persists 
qualitatively with a modulated Eckhaus line, while for large s, there is a qualitative 
difference to the unperturbed case that reflects the new network topology, i.e. the 
existence of a new loop. 

The article is structured as follows. In section 2 we discuss the stability of the 
zero solution (2), its spectrum and its bifurcations. In section 3 we find asymptotic 
expressions for the emerging periodic solutions in the case of small perturbation 
size s. We also reduce the case of asymptotically large s to the analysis of an 
inhomogeneous ring. Section 4 deals with the stability analysis of the periodic 
solutions and the results are compared to numerical simulations. Finally, we discuss 
our findings and give an outlook on possible extensions and applications of the 
results in section 5. 

2 Stability of the zero solution 

To study the stability of system (1) it is convenient to identify each variable Zj {t) G C 
with a two-dimensional real variable Zj (t) = {zj^i (t) , Zj ^2 {zj (t)) , A {zj G 

Then, system (1) is equivalent to the real system 



( 3 ) 


where 


is the representation of the multiplication with /i in 
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Figure 1: Coupling scheme of a unidirectional ring with a shortcut. 


2.1 Spectrum of the equilibrium 

We linearize system (3) in Zi = .. = zjv = 0 and obtain the variational equation 

it) = M^5zj + 5zj+i it) , j = 1,...,N - 1, 

(t) = M^Szn + Szi (t) + s5zi (t ), 

which can be written as 

j^6Z (t) = [Hat 0 + G, 0 Id2] (50 (t ), (4) 

where 6Z = (dzi, ...,6z]sf)~^, Mat G is the iV-dimensional identity matrix, 




"01 O' 

0 1 

1 0 s 0 


( 5 ) 


is the coupling matrix of the perturbed ring, A® B denotes the tensor product of 
the two matrices A and B. Equation (4) is a simple example of a system, which 
is treatable by a master stability function (MSF) approach [24], In our case the 
MSF M : C —)■ M simply reads M (A) = (U + (A), where A is an eigenvalue of the 

coupling matrix Gs- Indeed, the spectrum of (4) is 


a (Mat ® M^ + Gs® Id2) = a (M^) + a (GJ = {/i, /i} + (T (G^) . (6) 

Taking the real part, we obtain M (A) = (U + 3? (A). The spectrum of the coupling 
matrix Gs equals the set of solutions of the characteristic equation 

Xa. (A) = A" - sA'-‘ -1=0. (7) 


3 







Figure 2; Spectra of the coupling matrix Gs (5) for = 20 oscillators, a shortcut 
at node i = 6 and for different coupling strengths s: a) s = 0.1, b) s = 0.6, c) s = 1, 
d) s = 5. 


For s = 0 , cr (Go) consists of the N roots of unity 

• 27rfc _ 

7iv,fc = e ^ , k = 0, ...,iV - 1. 

For small s 7 ^ 0, the roots A^, k = 1 ,..., N, of eq. (7) are given (asymptotically) as 

Afc (s) = 7Ar,fc + ^ ) ( 8 ) 

as one can readily compute by applying the implicit function theorem to (7) with 
base points (sq, Aq) = ( 0 , 7 Ar,fc). Hence, the spectrum of Gs is a weak modulation 
of the spectrum of the circulant matrix Go [see Fig. 2(a), (b)]. For large s, a (G^) 
can be computed in a similar manner. In Appendix A we show that in this case the 
spectrum of Gg splits into two parts: there are ^ — 1 roots 

Ai,fc (s) ~ A: = 0, ...,£- 2, (9) 

located close to an inner circle of amplitude ~ and N — I + 1 roots 

A 2 ,fe (s) ^ fc = 0,..., iV - (10) 

which are close to an outer circle of amplitude ~ j^gge pig_ 2 ]. 

2.2 Bifurcations of the equilibrium 

From the formula ( 6 ) for the spectrum of the equilibrium Zi = .. = Zjv = 0, it 
follows that it is asymptotically stable iff 

3? (A) < —a, for all A G a (G^). (11) 

When the parameter a is increased starting from a value which satisfies (11), a 
bifurcation takes place whenever for some A G a (G^): 

a = -gfJ(A). (12) 
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Note that there is always one purely real eigenvalue Ai (s) = 1 + ;^ + (s^) of Gg 

which has maximal real part among all eigenvalues of Gg- Indeed, for any solution 
A = > 0 , of (7), we have 

1 = - sg^-^ = xg. (p) + 1- 

Hence, xGs (p) < 0, and lima,^.oo XGs (^) = +oo implies that there exists a real 
solution ^0 ^ ^ such that XGs (Po) = 0- Therefore, for small s > 0 the equilibrium 
switches stability at 



Since we assume /3 7 ^ 0, this bifurcation is a Hopf bifurcation and the emerging 
periodic orbit has frequency (3 at onset. As for s = 0 , the bifurcation is supercritical 
for small s > 0 , because the cubic term of the corresponding normal form depends 
continuously on s. Therefore, a branch of stable periodic solutions emerges and exists 
for a > «! (s). A further increase of a leads to a sequence of Hopf bifurcations which 
give rise to A^ —1 branches of periodic solutions. As in the case s = 0 , all bifurcations 
are supercritical if s > 0 is small (by continuity). The same is true if s > 0 is 
large, as we show in Appendix B. The latter A^ — 1 periodic solutions are initially 
unstable, inheriting instability from the steady state. The initial frequency cu of the 
emerging periodic solution equals the imaginary part of the crossing eigenvalue, that 
is cu = /3 + A (A) for the corresponding A G Gg. For the unperturbed ring (s = 0) the 
hrst L(A^ — l)/2j ([xj := max{? 7 , G Id : n < x}) branches stabilize when they cross 
the Eckhaus stabilization line [14] 



where |2’|^ is the (constant) amplitude of a periodic solution Z {t) = (zi (f),..., zn 
This observation is in striking analogy with the well known Eckhaus destabilization 
in diffusive systems [22, 23]. Remarkably, it is also found in this unidirectional sys¬ 
tem which does not extend to a spatially extended system in a natural manner. In 
section 4 we investigate how the added shortcut changes this scenario. 

3 Emergent periodic orbits 

Let Z (t) be a periodic solution of (1), which emerges from a Hopf bifurcation at 

a (s) = —3? (A (s)), (14) 

and which is associated to the eigenvalue A (s) of the coupling matrix Gg, see ( 8 )- 
(10). Because of the F^-symmetry of the system, we employ the ansatz 

Z (f) = (e, s), (15) 
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(a) k = l (h) k = 3 (c) s = 0.2 




Figure 3: Moduli of the components of the initial profiles V (0, s) = b (s) of emerging 
periodic solutions [see (17) and (19)] for different wave numbers k (i.e., different 
eigenvalues (s) of Gg) and coupling strengths s as indicated in the hgure. For all 
panels: N = 100 and i = 26. 


where 

e := a — a (s) > 0 (16) 

is the parameter distance from the bifurcation point, 

V (e, s) = (vi (e, s), ..., vjv (e, s))^ E C^, (17) 

is the profile vector and the frequency u (e, s) of Z [t) is 

a; (e, s) =/3 + 3^ (A (s)) + C> (e). (18) 

The emerging periodic orbit is e-close to the complex plane spanned by the eigen¬ 
vector b (s) of Gg corresponding to A (s) and tangential at the bifurcation point itself 
[25]. This means, V (0, s) = b{s) with 

b{s) = a (s) • (l, A (s), A^ (s),..., A^"^ (s))^ , (19) 

where one may assume a (s) G M due to the S'^-symmetry of the system. Figure 3 
shows several examples of the prohle shapes for different s and A. For |A| > 1 the 
emerging solutions become stronger localized at the iV-th node zn with increasing 
s since it scales as zn ~ For |A| < 1, the localization takes place at zi {t) for 

the same reason. 

3.1 Case I: small perturbation 

In this section we consider s to be small. Our aim is to determine a formal asymptotic 
expansion for the frequencies (18) and the profiles (17) of the periodic solutions and 
to derive evaluable approximate conditions for their stability. In particular, we are 
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interested in the deformation of the Eckhaus stabilization line (13). In order to 
obtain asymptotic expressions for the prohles in case that s,e > 0, we linearize 
the vector held in e = s = 0. For a periodic solution (15), we introduce rescaled, 
rotating coordinates Uj (t) E C, j = 1,N according to 



Zj = ^fee^^^^VjUj, 

(20) 

with u = u) ( 

e, s) and Vj = vj {e, s). Then (1) becomes 



iij = [a + i {13 — u) — e \vjUj\'^') Uj + 

'^j 

(21) 


un = {a + i {(3 — oj) — e IvNUNlf um + —Ui + s—u^, 

' Vn Vn 

(22) 

In rotating coordinates, the equilibrium solution 



Uj = 1, j = 1, 

(23) 

corresponds to the periodic solution (15) of (1) and the stability of (15) and (23) is 
the same. In Appendix D we show that for each eigenvalue (s) = 'yN,k + flN k + 
O {sf of Gg, the corresponding branch of periodic solutions has frequencies 


iUk{£,s) — /3 + A (7Ar,fc) + —(7v,fc) + C (( £ + s )^) , 

(24) 

and prohles 

Vj {e, s) - 7jv^^ ^ N ^ ^ ^ ^ ) • 

(25) 


3.2 Case II: large perturbation 

In this section we consider s to be large. To treat (1) as a weakly perturbed system 
we perform a change of variables 

yj (i) = > (26) 

with a small parameter which is the inverse radius of the outer spectral 

circle of eigenvalues (10). This transformation normalizes the emerging prohles (19) 
corresponding to the eigenvalues (10). The transformed variables (26) satisfy 

yj W = lyj (t)n Vj (t) + <^‘^^~^yj {t ), j = 1, ...,N - 1, 

i/N (t) = - \yN (t) n yN (t) + (t). (27) 

Since {‘^^^~^yi) and is small, we consider system (27) as a small 

perturbation of 

yj (t) = \yj (t)l^) yj {t) + {t) , j = 1, ...,N - 1, 

ijN if) = - \yN (i) n yN f) + f) ■ (28) 
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Although we cannot show that results for the persistence of hyperbolic invariant 
manifolds [26] apply and assure that (28) possesses the same hyperbolic invariant 
manifolds as does (27), the truncated system (28) is a natural approximation to (27). 
Note that in (28) the components yi, do not couple back to the rest of the 

system since the weak link from yi to yj^ was taken out. Therefore, the dynamics 
of the subsystem yi, ..., y^ is independent and, apart from the zero solution, acts 
as a periodic force on the attached subsystem yi, ..., yi-i. In original variables (28) 
reads 


Zj{t) = {y-\zj{t)f) Zj + , j = -1, 

zn (t) = {y-\zN {t)f) zn + sze (t) . (29) 

The linearization of (29) at the zero solution has eigenvalues 

fi + i', and ft + u 

where z/ is a root of the characteristic equation 

_ s) z/'-i = 0 

of the reduced coupling matrix Hg, which is obtained by erasing the link from zi to 
Zjsf from Gg. It has eigenvalues 

1 ^ = 0 and u = k = 1,N — i + 1. (30) 

Periodic orbits which correspond to the algebraically {i — l)-fold eigenvalue = 0 
are localized on the nodes Zi, ..., zi_i. 


4 Stability of the periodic orbits 

4.1 Case I: small perturbation 

The variational equation for system (21)-(22) at the stationary solution (23) is 


u 


— (cr E i (/3 — cu'j 

iiN = [a + e + i — cj) 


o I |2\ , '^i+l 

- 2e \vj\ ) Uj + -—Uj+i, 

I i 2 \ "^1 

— 2e uat ) un H- ui + s — ui. 

Vn Vn 


We transform the system into 
(24) and (25) to obtain 


p2Af 


X 


= ‘^{uj))'^) and insert the expansions 


Xj — — [-^A A Xj + M^Xj^i + (D (d^l + |s|) ^ 

ijv = -[M-^ + 2edii]xN + M-^Xi + sM)^i^[xi-Xi] + 0{{\e\ + \s\f) (31) 



with Ao = 7Ar,fc, A (s) = Aq + ;§^Ao, the matrix representation M : C —)■ [as in 

(3)] and 5mn = {^jm^kn)j fc IS the matrix which has the entry 1 at position (j, k) and 
zeros everywhere else. We drop the higher order terms in (31) and write the system 
in the form 


X — A. (£■, s) X — —Id/v" ® 2eSii) + Gq (E) + (^vi — ^Ni) ^ ^■ 

(32) 

Clearly, an MSF approach as in section 2 is not feasible anymore. However, we have 
reduced the dynamical problem to an algebraic one. The eigenvalues of system (32) 
can be computed by standard numerical procedures to determine approximately 
the stability of the periodic solutions in the vicinity of a bifurcation point. The 
eigenvalues of (32) approximate the eigenvalues of the exact system (21)-(22) at the 
steady state (23) up to hrst order in e and s. Anyway, this hrst order approximation 
leads to good predictions even for moderate values of the parameters, in particular 
for e [see Fig. 4]. 


Resonances 

An important observation for the perturbed system is that in comparison to the 
unperturbed case s = 0 , for small s > 0 the point of stabilization of a periodic 
solution may be altered or not. It remains nearly the same if the corresponding 
eigenvalue A^ = + O (s) fulhlls 

arg ( 7 v,fc) ~ arg ( 7 ^, (33) 

or equivalently, G Z. This corresponds to a situation where both components 

of the input Z\ {t) + sz^ {t) to node z^ {t) possess approximately the same phase. At 
the same time it is a condition for maximizing the modulus |Afc| ~ |7v,fc + k\ 
and for 77 V,a: to span an eigenmode of both, the unperturbed cycle of length N 
and of a cycle which has the same length N — i — 1 as the newly created cycle 
i ^ N ^ {N — 1) ^ ^ i, i.e. to be an At-th and an {N — i + l)-th root of 

unity. On the contrary, the antiphase condition 

arg ( 7 v,fc) ~ arg ( 7 ^ + tt. (34) 

causes the point of stabilization to grow rapidly with increasing s, i.e., it destabilizes 
the corresponding periodic orbit. The more precisely the equality (34) holds, the 
more pronounced is the destabilizing effect of the additional link [see hg. 4]. 


4.2 Case II: Large s 

For large s both systems, the original (1) and the truncated (29), admit two types 
of periodic solutions emerging in bifurcations corresponding to eigenvalues of scale 
|A| 0 or |A| ~ respectively [see eq. (9), (10), and (30)]. In case of system 
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Figure 4: Bifurcation diagrams for different strengths s of the shortcut in a ring of 
N = 100 oscillators with a shortcut from node £ = 26 to node N = 100; /3 = 2.5. 
Dashed lines indicate unstable periodic orbits and solid lines stable periodic orbits. 
In black the stabilization line for s = 0 is shown. The frequency of the solutions is 
plotted against the bifurcation parameter a. The perturbation strength s increases 
from left to right: a), d) s = 0.05; b), e) s = 0.1; c), f) s = 0.2. The upper panels a), 
b), c) show the approximated diagram obtained from (32). The lower panels d), e), 
f) show results of numerical bifurcation analysis of the full system for comparison. 
These calculations were carried out with the program AUTO [27]. 


(1) all these periodic solutions emerge in form of rotating eigenvectors (15)-(19) of 
the coupling matrix Gg which leads to a locally pronounced activity in zi (t) for a 
corresponding eigenvalue |A| ~ 0, and in zn (t) for |A| ~ g^^-e+i rjpg 

same picture applies for bifurcations of (29) corresponding to the simple eigenvalues 
(|z/| of Hg, while the i — 1-fold bifurcation at a = 0 corresponding to 

iz = 0 simultaneously creates several solutions which are completely localized in the 
attached subsystem {zi, ...,z^-i) where only the oscillators zi, ..., Zk, k = 1, —1, 

are active and all others silent. For each of these solutions the first active element 
Zk {t) = is located on the limit cycle of an isolated Stuart-Landau oscillator 

and all other Zj, j < /c, lock either in phase or antiphase to their input signal 
Zj+i- However, these solutions can never stabilize, since the zero solution of the 
subsystem zi, ..., zn is unstable after the hrst bifurcation corresponding to |z/| 
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Therefore, it suffices to study the inhomogeneous ring 

Zj {t) = (/i - \zj {t)f) Zj + Zj+i {t), j = 1 ,n - 1, 

Zn{t) = {jj,- \Zn{t)f) Zn + SZi{t) , n = N - i +1, (35) 

to understand the possibly stable dynamics of (29). To approximate the stability of 
the emerging periodic solutions one can proceed as for the case of small s: write the 
system in scaled rotating coordinates (20), linearize around the equilibrium solution 
(23), expand the variational equations in powers of e and truncate terms of order 
higher than O (e).We obtain the following approximate variational equation [see 
Appendix E]: 







3 

0 


0 

1 





1 0 
0 1 




(36) 


I 2 2(_j' — 1) — 

where |u°(s)| = ns ^ Although the loss of symmetry prevents us from 

applying an MSF approach, (36) enables us to approximate the Floquet exponents 
by solving numerically the characteristic equation of (36). 


5 Discussion 

We have investigated how the introduction of a shortcut alters the dynamics in 
a unidirectional ring of Stuart-Landau oscillators. In absence of a shortcut, the 
system exhibits a bifurcation scenario similar to the Eckhaus instability observed in 
dissipative media. For small shortcut strengths s we have found that the Eckhaus 
stabilization line is modulated in the following manner: The destabilizing impact 
on periodic solutions is stronger for non-resonant modes than for resonant ones. 
The latter correspond to wavenumbers that are compatible with the lengths of both 
cycles which exist in the perturbed system, i.e. for the cycle of length N — i 1 
which contains the shortcut, as well as for the full cycle of length N. In contrary 
to the non-resonant solutions, the stabilization of the resonant periodic solutions 
occurs for similar parameter values as in the case without shortcut. As a result, one 
can control the destabilization of a specihc set of wavenumbers via the link position 
i and its strength s. In the case of a large shortcut strength s we have provided 
an argument that the cycle of length N — i + 1 dominates the dynamics and stable 
solutions can be treated as solutions of a single unidirectional inhomogeneous ring 
which has coupling strength s at only one link. 

Further investigations will be dedicated to how small perturbations may be used 
to select solutions with a specihc wavenumber by adding a corresponding set of res¬ 
onant links. More generally, our Endings may help to understand how perturbations 
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0) 


Figure 5: Bifurcation diagram for intermediate strength s = 5 of the shortcut. Fixed 
parameters are N = 100, i = 26, (3 = 2.5. Branches of unstable periodic orbits are 
indicated by dashed lines and stable branches by solid lines. Along the black line 
(13) the stabilization takes place for s = 0. The frequency of the solutions is plotted 
against the bifurcation parameter a. The calculations were carried out with the 
program AUTO [27]. 


with a more complicated structure, consisting of several shortcuts, can influence the 
dynamics of a unidirectional ring. Moreover, our observations might even help to 
locate an unkown shortcut when one is only allowed to vary a bifurcation parameter 
and observe the dynamics, since the modulated Eckhaus line can be used to identify 
the shortcut. A strong shortcut can be used in arbitrary networks in order to local¬ 
ize the activity on the cycles in which they are contained and amplify the activity 
in particular on their targets. 

Another point which was not investigated here but deserves a closer study is 
the development of prohles far from the bifurcation point. The simplest, most 
important phenomenon is that, for increasing «, the exponential tails of the prohles 
develop into plateaus, where the prohle amplitude is locally constant as a function 
of the component index. By taking the limit N ^ oo one can argue that solutions 
may consist of several, sharply separated plateaus where each plateau can possess 
a different wavenumber. This complies with numerical observations, although all 
observed solutions with more than one plateau lie on the unstable branches which 
correspond to the inner spectral circle for large s. For sufficiently large values of s, 
those branches begin to curl with increasing a and seem to be unstable for arbitrary 
large values of a [see Fig. 5]. 
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A Adjacency spectrum for large s 

To apply the implicit function theorem and continue roots of the characteristic 
equation (7) for the case of large s, we define 

= tX^ - (37) 

Then F (A, -i, i) =0 is equivalent to (7). We now apply the implicit function 

theorem twice to find the two distinct families (9) and (10) of small and large 
solutions of (7). Let us first compute 

dxF{X,T,^) = tNX^-^ -{i-l)X^-\ 

drF{X,T,^) = X^, d^F{X,T,^) = -1. (38) 


Consider the eqnation 


F (A, 0,^9) = -X^-^ -^ = 0. 


(39) 


It possesses ^ — 1 solutions 


Ai,fc (0,^9) = 

A; = 0, — 2 . Assuming (we show that below) that for (r, ■d) 7 ^ 0 one can extend 

these solutions to smooth functions (r, d) Ai^^ (r, ■d) which solve F (Ai^^ (r, d) , r, d) 
0. These can be expanded in r = 0 as 

Ai,fc (r, d) = Ai,fc ( 0 , d) - dxF (Ai,fc ( 0 , -d), 0 , h*)”^ drF (Ai,fc ( 0 , -d) ,0 ,'&)t + O (r^) 

= Ai,, (0, ^) + ((£ - 1) X{J (0, ^)) A(^, (0, ^)t + 0 (r^) 

= Ai,fc( 0 ,'d)H-^^-r 

and, for d = t, 


- 1 


+ O (r^) 


Ai,fc (r, r) = Ai,fc ( 0 ,r) + 


Ai!,-^+^(0,r) 

£-1 


+ O (r^) 


(40) 


This gives a family of £ — 1 solntions of (7) situated near a small circle of radius 
1 

~ (l/s)^-L For r 7 ^ 0 , the equation 


F(A,r,0) =rA^-A^-^ = 0 


(41) 


has an (£ — l)-fold root at A = 0 which corresponds to the solutions of (39) and 
N — i + 1 roots ^ 

A 2 ,fc (r, 0) = |r|"^-^+i 7 w-£+i,fc (42) 

A; = 0 ,As before, we obtain an asymptotic representation 
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A 2 ,fc(r,r) = A 2 ,fc (r, 0 ) - 5aF (A 2 ,fc (r, 0 ), r, 0 ) ^ (A 2 ,fc (r, 0 ), r, 0 ) r + C> (r^) 

= A2,fc(r,0) + XG, (A2,fc(r,0))"V + C>(r^) (43) 

This is the family of solutions which lie near a larger circle of radius ~ . 

Now let us show that the interval of existence of the implicit functions r 
Ai,fc ('7‘, "do) and -d i—)■ A 2 ,fc ('To^'d), respectively, indeed contain the points r = -do and 
■d = To, respectively. Let us denote the implicit function in question simply A (r, d). 
If the implicit function theorem fails to provide an extension of A (r, d) in some point 
> 0 , we must have 

dxF (A*, n, d*) = - 1) Af-" = 0 (44) 


where A* = A (r*,^*). Since F (0,r, d) = 0 is equivalent to d = 0, we may assume 
A* 7 ^ 0. Thus, (44) is equivalent to 


^N-e+i 


i-1 

1 ^' 


(45) 


Furthermore, from (37) 
(44) to obtain 


0 we obtain r* = A^ ^ ^ ^ which we insert into 




N 

N-i + 1 


■d*. 


(46) 


From (45) and (46), we obtain 


r* = F (d*) 


i-l 

N 




N-^+1 




N-£+l 

£-1 


Since F : M>o M>o is monotonic, we have 


dxF (A (r,d*) ,r,d*) 7 ^ 0 

for all r < r* = F (d*). Since F (r) —)■ cxo, for r \ 0, there exists Tq > 0 such that 
r < F (r) and A (r, r) is dehned uniquely, for 0 < r < To. 


B Supercriticality of the Hopf bifurcations 

To show that the bifurcations at a = —3? (A) of system (3) are supercritical for 
sufficiently large s > 0, we use the projection method for center manifolds [25]. In 
the following, we write the vector field of (3) as 

fa{z) = Az + lc {z,z,z), 
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where A = Mat (8) + Gg® Id 2 is the linearization of /q, at z = 0 and the trilinear 

function C contains all cubic terms. Furthermore, let v G be an eigenvector of 
A corresponding to the eigenvalue /i +A, A G o'(Gs). (The case of an eigenvalue fi + X 
can be treated analogously.) Let w G be the normalized adjoint eigenvector 
corresponding to v, i.e. = (/f + X)w and {w, v) = 1. We have 


V = (l,A,...,A^ ^)^(g) (47) 

w = (48) 


with K = 2A^ — 1) + {N — ^ + 1) X^). A Hopf bifurcation at a is supercritical 

for a negative and subcritical for a positive hrst Lyapunov coefficient 

( 0 ) = > ( 49 ) 

/Wq 

where wq = /9 + O' (A) 7 ^ 0. Writing C = (Gi^i, Gi_ 2 , • • •, Cn, 2 ) we obtain 

Gj -1 {v, v,v) = -6vj^i - 2 {2vj^i \vj^ 2 \^ + t'|, 2 Wi) = iCj ^2 (v, v,v) , (50) 


for 1 < j < A^. Using this, we get 


{w,C{v,v,v)) 



(£ - 1) + (A^ - £ + 1) A^ 


(51) 


For large s > 0 we distinguish two cases.For the case |A| ~ gnd that 

'Si{{w,C {v,v,v))) ^ — 00 , as s ^ cxo, 

and for the case |A| ~ 

^{{w,C {v,v,v))) 0, as s — )> 00 . 

This means that for sufficiently large s, we have li (0) < 0 for all eigenvalues. Thus, 
all bifurcations are supercritical. 


C Supercriticality for an inhomogeneous ring 


Consider the case 1=1 and s arbitrary, i.e. a ring with inhomogeneous coupling 
strengths. Here, using (51) and the characteristic equation, X^ = 1 + s, we get 


(w,C(v,v,v)) 


8(l-(l + .sf) 

jV (1 - |1 + s|’''^) 


(52) 


Thus, in the case of the inhomogeneous ring the bifurcation is supercritical for 
arbitrary s. 
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D Expansion of the solution profiles for small per¬ 
turbations 

The linearization of (21)-(22) at (23) is 

0 = [a-\-i {13 — u) — (53) 

'^3 

0 = [a + i{l3 - u) - e\vNf‘) + — + s—, 

Vn Vn 

From these equations we will now find expressions for the first terms of the Taylor 
expansions of the unknown functions 

(^; s) = CJoo + ^<^10 + ^ ((kl + kl)^) 

and 

Vj {e,s) = + evf + svf + O ((k| + |s|)^) , j = 1, 

Terms at order O {1) 

Considering the terms in s = e = 0 in (53) yields the circular equations (let := 


yOO 

0 = + i(/3 — Woo) H—(54) 

yOO 

with the shorthand (0) = — cos {27rk/N). This leads to 

(-a° -i{/3 - woo))^ = 1 

which contains no new information, since it only determines = —3? (Aq) and 
1^00 = k + ^ (Aq) with an A^-th root of unity 

, , ^ 2Tvk 

Ao = Afc (0) = e* ^ . 

For the profile, (54) yields 

= Xiro, (55) 

with a hitherto unknown scale factor Vq := n™. Without loss of generality, one may 
choose Vi {s,s) G M>o, because of the phase shift invariance of (15). In particular, 
we then have tq G M+. 


16 




Terms at order O (e) 

At first order in e we obtain another set of circnlar eqnations (let again 


0 = 1- iuio - ’^0 + 
which leads us to the following recursion 


A, 


0 I ^j+i 


,10 


^0 V Xi 




V 


10 

'1+1 




10 


(56) 


Dehning aj = -~=rr and A = (rQ — (1 — iuio)) equation (56) can be written as 


Ao' 


Oj+i = A + ttj with the solution Oj+i = jA + ai. Because of the circularity otv+i = ni 
we then have oi = NA + ai which determines A = 0. Therefore, Tq = (1 — iujio) and 
hnally 

wio = 0 and ro = 1 (57) 

That means at hrst order the frequency of the oscillations does not depend on e. 
For the prohles it follows (as at order O {!)) 


’15 


= AJrio, 


(58) 


with rio := Ui G M. 


Terms at O (s 


Due to the fact that s perturbs the network symmetry, the equations at order O (s) 
are more complex than those at order O (e). Again, the linear terms in s of (53) 
give a recursive formula 

,01 


^1+1 ^ 

K 

= 


— (zwoi - a ) + 

^0 An 


Ar 


{iujQi - - Aq ^ + 


V 


01 

N 


, V-1 ’ 


(59) 


'0 


yUL 

where = —^3fJ(Afc(s)) |s=o- Foi’ j < N we have {ioooi — a^) + 


Aq 


Inserting the resulting expression for in the second equation gives = 
^ {iuoi — a^) — Ag”^ + Therefore, N (iwoi — «^) = Ag or equivalently 


^01 = (Ag) and 3? (A^ . 


N 


The perturbations Vj]_i of the prohle are then determined as 


A. 


^ + ’’01 ) K- 


(60) 


(61) 


up to a scaling rgi := G M as above. 
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Higher order terms 

To determine the amplitudes rio and roi of the perturbations and j = 
1, we need to calculate the second order terms O (e^), (9 (s^), and O {e ■ s) of 

(53) due to the nonlinear term. We omit this here and just give the resulting values 

rio = 0 and rgi = 0. (62) 

The vanishing rio means that at hrst order the prohles do not depend on e. 


E Expansion of the solution profiles for the inho¬ 
mogeneous ring 


A periodic solution of (35) corresponds to a hxed point in co-rotating coordinates 
(20). The stability of this hxed point is governed by its variational equations 

Xj = ^Md+idi-uj) - £ (s) \vjf‘ 

Xn oj) £ (®) |^n| 

Higher order terms of u = oj (e, s) and Vj = Vj {e, s) can be determined by the 
equations 

0 = {a + i {P — oj)) — e \vjf -{ —(65) 

'^3 

0 = {a + i{l3 — uj)) — e\vnf‘+ s—, (66) 

'^n 


3 0 
0 1 

'3 0 
0 1 


Xj + M^jjj^Xj+i, 

Xn + sMn^Xi, 


(63) 

(64) 


which are obtained from inserting the solution ansatz (15) into (35). Solving for real 
and imaginary parts yields the conditions 



(67) 

( 68 ) 


Note that we expand the unknown functions only in e, keeping s arbitrary. Using 
(67) and (68) the variational equations (63)-(64) write 


Xi 


Xn 


M.j^ + 2£ (^ J JJ 

sM^+2e\v„f (^l 0 


Xj + M^^Xj+i, 

'^3 

Xn + sM^Xi. 


(69) 

(70) 
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In the rest of this section, we hnd expressions for (0, s) and ^ 
which can be inserted into (69)-(70) to obtain (36). Let 



|e=0 


uj [e, s) — cuq (s) + EOJi (s) + O (s'^) 


and 

Vj (e, s) = (s) + EVj (s) + O (e^) , j = 1, N. 

For e = 0, equations (65)-(66) yield 

Vj+i = {i {uo -13)- ao) = ... = {i{uo -13)- aoY n?, 

Vi = ^ {i {uq -/3) - ao) v°, 

where Oq = —3? (A^ (s)) = —cos {27rk/n) is the critical value at which the periodic 
solution emerges, ujq = (3 + (A^) is the initial frequency, and the initial prohle is 

given by Vj = 

For terms of order O (e) one obtains 


cui = 0, 




- 1 


= n- 


— 1 


and v] = {\v^_i\^ - l) + A^ 

1=0 




To approximate (69)-(70), we calculate the hrst order terms of the quotients Vj+i/vj, 
which gives at 0{e^) 


0 /n,/c? 

yO 


and at 0{e) 


^i+i 


1 



- 1 . 
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